One-dimensional Disordered Density Waves and Superfluids: 
The Role of Quantum Phase Slips and Thermal Fluctuations 
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The low temperature phase diagram of ID disordered quantum systems like charge or spin den- 
sity waves, superfluids and related systems is considered by a full finite T renormalization group 
approach, presented here for the first time. At zero temperature the consideration of quantum 
phase slips leads to a new scenario for the unpinning (delocalization) transition. At finite T a rich 
cross-over diagram is found which reflects the zero temperature quantum critical behavior. 
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The collective behavior of condensed modulated struc- 
tures like charge or spin density waves (CDWs/SDWs) 
§, |L flux line lattices (FLLs) § § and Wi gner crys- 
tals P] in random environments has been the subject of 
detailed investigations since the early 70th. These are 
motivated by the drastic influence of disorder: without 
pinning CDWs would be ideal superconductors whereas 
type-II superconductors would show a finite resistivity. 
In three dimensional systems the low temperature phase 
of these structures is determined by a zero temperature 
disorder fixed point resulting in quasi-long range order 
and a glassy dynamics [|[ || . In two dimensions this fixed 
point is extended to a fixed line which terminates at the 
glass transition temperature |6j. In the low tempera- 
ture phase correlations decay slightly faster than a power 
law and the linear resistivity vanishes (for a recent review 
see Q). In one dimension the situation is different: the 
glass temperature is shifted to T — 0. Nevertheless there 
remains a residual trace of the disorder which is reflected 
in the low temperature behavior of spatial correlations 
and the dynamics (7j. 

Clearly, at low temperatures also quantum fluctuations 
have to be taken into account. Disorder and quantum 
fluctuations in ID CDWs at zero temperature have been 
considered previously (see e.g. [^], ^]) and an unpinning 
transition as a function of the interaction strength was 
found. Finite temperature effects were partially incorpo- 
rated by truncating the renormalization group (RG) flow 
at the de Broglie wave length of the phason excitations ||. 
However, for a complete study of the thermal to quantum 
cross-over quantum and thermal fluctuations have to be 
considered on an equal footing Jl(| , which will be the first 
goal of this paper. The second goal will be the consid- 
eration of quantum phase slips which trigger their own 
quantum phase transition to the disordered phase and 
influence also the unpinning transition. Experimentally, 
ID behavior can be seen in real materials e.g. in whiskers 
with a transverse extension smaller than the correlation 
length or in chain like crystals with weak interchain cou- 
pling. In the latter case there is a large cross-over length 
scale up to which ID behavior can be observed R ||. 



The results obtained for the CDWs/SDWs have a number 
of further applications on disordered quantum systems: 
they relate e.g. to the localization transition of Luttinger 
liquids || [|, tunnel junction chains |pd| , superfluids |fl2| , 
Josephson coupled chains of these systems if the coupling 
is treated in mean-field theory || . In most parts of this 
paper, we will use the terminology of CDWs. 

Well below the mean field condensation temperature 
Tmf of the CDW, the electron density p(x) can be writ- 
ten in the form ffl 



p(x) = p (l + Q 1 d x <p) + pi cos[p(t / 3 + Qx)} 



(1) 



where Q = 2k p, kp is the Fermi-momentum, po the 
mean electron density and p\ is proportional to the am- 
plitude of the complex (mean field) order parameter 
Ae l,p ~ (6q + &1q). b\ , bk, denotes the phonon creation 
and annihilation operator, respectively. tp(x) is a slowly 
varying phase variable. Neglecting fluctuations in A, the 
Hamiltonian of the CDW is given by 



H = 
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+ITcos(^ J dyP{y)) \dx {2\ 



f&r* denotes 



where [P(x),<p(x')] = j5(x - x'). c 
the elastic constant, Vp the Fermi-velocity, v the effec- 
tive velocity of the phason excitations and f(T) the con- 
densate density @. Note that f(T) and A(T) vanish 
at Tmf whereas v remains finite. The third term re- 
sults from the effect of impurities of random potential 
strenght Ui = ±Ui mp and position xi and includes a 
forward and a backward scattering term proportional to 
Po and p±, respectively. We will assume, that the mean 
impurity distance limp is large compared with the wave 
length of the CDW and that the disorder is weak, i.e. 
1 *C hmpQ <C c/(Ui mp pi). In this case the Fukuyama- 
Lee length Lfl = {c/(Up 2 )) 2 ' 3 is large compared with 
the impurity distance, here U = Ui mp p\/ \ 
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fourth term in (@) describes the influence of quantum 
phase slips by tp = ±qit and will be further discussed be- 
low E3|. The model (0) includes the four dimensionless 
parameters t = T/irAc, K = hv/irc, u 2 — U 2 /A 3 irc 2 and 
w = W/ircA 2 , which measure the strength of the thermal, 
quantum, and disorder fluctuations and the probability 
of phase slips, respectively. A = ir/a is a momentum 
cut-off. Although for CDWs and SDWs if-values of the 
order 10 _1 and 1, respectively, have been discussed at 
T = [^3|, the expressions relating K and t to the mi- 
croscopic theory lead to the conclusion that both diverge 
by approaching Tmf whereas the ratio K/t remains fi- 
nite. The classical region of the model is given by K <C t 
which can be rewritten as the condition, that the ther- 
mal de Broglie wave length At = K/ (tk) of the phason 
excitations is small compared with the lattice spacing a. 

In order to determine the phase diagram we adopt a 
standard Wilson-type renormalization group calculation, 
which starts from a path integral formulation of the par- 
tition function corresponding to Hamiltonian (||) with 
u,w <C 1. We begin with the renormalization of the dis- 
order term and put w — for the moment. The system is 
transformed into a translationally invariant problem us- 
ing the replica-trick. Going over to dimensionless spatial 
and imaginary time variables, Aa; — > x and Avt — > r, the 
replicated action is given by (a — (upoA/ piQ) 2 ) 



1 
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dr' 



u 2 cos 



pUp a {x,T) - <Pp(x,T') 



+ ad x cp a (x,T)d x (pp(x,T') 



(3) 



Integrating over the high momentum modes of (p(x, r) in 
a momentum shell of infinitesimal width 27r/& < \q\ < 2n 
but arbitrary frequencies and rescaling x — > x = x/b, 
t — > 7~' = r/b we obtain the following renormalization 
group flow equations (up to one loop): 



dK 
~dT 

dl 



--p\ 2 KB a (p 2 K,-)coth-, 



p 2 K , K 
—— coth — 
2 2t 



dt _ 



(4) 
(5) 
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Bi{v,y) = dr 







dxgi (r, x) cosh (y — t) (cosh y) 1 



1 + (f) 2 (^sh f- COS ^) 



v/A ' 



(6) 



where I = In 6 and go(r,x) = 6(x)r . Note that 

K\ 
It > 



B (p 2 K, f ) for K -> 0. 



The equation for the flow of a is more involved and 
will not be discussed here since it does not feed back 
into the other flow equations. Indeed, we can get rid 
of the forward scattering term by rewriting <p(x) — 
ip b (x) +(pf(x) with <ff(x) = f*dyc(y), (c(x)) = and 
(c(x)c(x')} — ^aS(x — x'). The phase correlation func- 
tion C{x,t) = (((p(x,T) - (ys(0, 0)) 2 ) = C(x,t) + C f (x) 
has therefore always a contribution C/(x) ~ |x|/^/ with 
ZJ 1 ~ cr(l = log I x\). Since all further remarks about 
phase correlations refer to Cb{x, r) we will drop the sub- 
script b. 

There is no renormalization of t (i.e. c) because of 
a statistical tilt symmetry [[l4|. The special case t = 
was previously considered in || (with p = \/2). The 
flow equation for K obtained in H] for w = t = devi- 
ates slightly from (^), which can be traced back to the 
different RG-procedures. The critical behavior is how- 
ever the same: there is a Kosterlitz-Thouless (KT) tran- 
sition EE) at K u between a disorder dominated pinned 
and a free unpinned phase which terminates in the fixed 
point K* = 6/p 2 . uq denotes the bare value of the disor- 
der and K u is given by u 2 = ^( ^5" - log ^f-) with 

rj = Bo(p 2 K*, 00). In the pinned phase the parameters 
K, u flow into the classical, strong disorder region: K — > 
0, u — * 00. The correlation function C(x, 0) ~ |x|/^ M 
increases linearly with Integration of the flow equa- 
tions gives for small initial disorder and K <C K u an 
effective correlation length £ u w A _1 (ALfl) ( ' 1 ~ k I Kv ^ 
at which u becomes of the order unity. Close to the 
transition line £ u shows KT behavior. For K > K u 
£ u diverges and C(x,t) ~ K(l = log \z\) log |z| where 
|z| = Var^+r 2 . Note that saturates on large scales 
at a value K c g(uo). 

For large values of u our flow equations break down, 
but we can find the asymptotic behavior in this phase 
by solving the initial model in theit strong pinning limit 
C/i mp 3> 1,K — exactly. A straightforward but 
somewhat clumsy calculation yields for the pair cor- 
relation function C(x,t) = ||(1 - ^^)\Qx\ where 
a = tt / (pQl i mp ) . The connection to the weak pinning 
model follows by choosing Zj mp sa ^fl- 

At finite temperatures thermal fluctuations destroy the 
quantum interference effects which lead to the pinning of 
the CDW at t = 0. The RG flow of u in the region 
K < K u first increases and then decreases. £ can be 
found approximately by integrating the flow equations 
until the maximum of u(l) and t(l)/(l + K(l)) is of or- 
der one. This can be done in full generality only nu- 
merically (see Fig. 1). It is however possible to discuss 
several special cases analytically. The zero temperature 
correlation length can still be observed as long as this 
is smaller than which rewrites for K not too close 

to K* as t < t K w Kt { u~ K,Ku) ~\ t u « (ALfl)- 1 . 
We call this domain the quantum disordered region. For 
K > K u the correlation length £ is given by At which 
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is larger than given by purely thermal fluctuations. For 
scales smaller than At, C(x,t) still increases as ~ log \z\ 
with a continuously varying coefficient K e g(uo). In this 
sense one observes quantum critical behavior in that re- 
gion, despite of the fact, that the correlation length is 
now finite for all values of K [ju^. In the classical dis- 
ordered region tx < t < t u the correlation length is 
roughly given by Lfl as follows from previous studies 

|?J . In the remaining region t u < t we adopt an alterna- 
tive method by mapping the (classical) one-dimensional 
problem onto the Burgers equation with noise ]l6| ]. In 
this case the RG-procedure applied to this equation be- 
comes trivial since there is only a contribution from a 
single momentum shell and one finds for the correlation 
length £" x sa f /(T)t(l + (2tt /p) 2 (t u /t) 3 )A. The phase 
diagram depicted in Fig. 1 is the result of the numerical 
integration of our flow equations and shows indeed the 
various cross-overs discussed before. 

So far the phase field was considered to be single val- 
ued. Taking into account also amplitude fluctuations of 
the order parameter the phase may change by multiples 
of 27r by orbiting (in space and imaginary time) a zero 
of the amplitude. Such vortices correspond to quantum 
phase slips described by the last term in (||) (with q = 2), 
which we discuss here under equilibrium conditions. This 
operator superposes two translations of (p by ±qir left 
from x, i.e. it changes coherently the phase by itpr in 
a macroscopic region. For vanishing disorder the model 
can be mapped on the sine-Gordon Hamiltonian for the 
0-ficld (with K replaced by K~ r ) by using the canoni- 
cal transformation P = — ~d x and —~d x <p = Id. To see 
the connection to space-time vortices one rewrites the ac- 
tion of interacting vortices as a classical 2D Coulomb gas 
which is subsequently mapped to the sine-Gordon model 
p7|. The initial value wq of w is proportional to the fu- 



jacity wo 



-S core 



l h of the space-time vortices which 



may be non-negligible close to Tmf , where the action 
Score ~ h/(irK) of the vortex core is small. Performing 
an analogous calculation as before (but with u = 0) the 
RG-flow equations read: 



dK 

Hi 
dt 
~dl 

dw 

~dJ 



-^B 2 f^,-)coth^, 
2 K 3 \K 2tJ 2i' 



K 



2 K 4 \K 2t 



coth • 



K 

Yt 



^— coth — 
4K 2t 



w, 



(7) 
(8) 
(9) 



where are given in (||) with g\ — 2r 2 cosa; and 

g 2 = (x 2 + t 2 )cosx. From <Mj - (0) we find, that for 
t = u = quantum phase slips become relevant (i.e. w 
grows) for K > K w with #* = q 2 /8 (q = 2 for CDWs). 
In this region vortices destroy the quasi long range or- 
der of the CDW, C{x,t) ~ \z\/£, w . The transition is of 
KT type with a correlation length £ w (w(log£ w ) ~ 1) 
diverging at K w + |L2| . At finite temperatures w first 



increases, but then decreases and flows into the region of 
large t and small w. Thus quantum phase slips become 
irrelevant at finite temperatures. This can be understood 
as follows: at finite t the ID quantum sine-Gordon model 
can be mapped on the Coulomb gas on a torus of perime- 
ter Kj t since periodic boundary conditions apply now in 
the T-direction. Whereas the entropy of two opposite 
charges increases for separation L ^> K/t as log(LK/t), 
their action increases linearly with L. Thus, the charges 
remain bound. The one-dimensional Coulomb gas has 
indeed only an insulating phase (L8| . 
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FIG. 1: The low temperature cross-over diagram of a one- 
dimensional CDW. The amount of disorder corresponds to a 
reduced temperature t u ~ 0.1. In the classical and quantum 
disordered region, respectively, essentially the t = behavior 
is seen. The straight line separating them corresponds to 
At ~ a. In the quantum critical region the correlation length 
is given by At- Pinning (localization) occurs only for t = 
0,K <K*. 




FIG. 2: T = phase diagram for a CDW with quantum phase 
slips. If qp < 4\/3 there is a single transition between a low- 
K pinned and a high-K unpinned phase. In both phases the 
correlation length is finite. If qp > 4\/3 these two phases are 
separated by a third phase in which phase slips are suppressed 
and C(x,t) ~ log \z\. Both transitions disappear at finite t. 

It is now interesting to consider the combined influence 
of disorder and phase slips. In doing this we write an 
approximate expression for the action of a single vortex 
in a region of linear extension L as 



Sv 



2)logL 



K ' 



(10) 



For very low K{< K U ,K W ) where u(L) w uoL 3 / 2 the 
disorder always favours vortices on the scale of the ef- 
fective Fukuyama-Lee length £ u . These vortices will be 
pinned in space by disorder. On the other hand, for 
very large values of K(> K U ,K W ) phase vortices are 
not influenced by disorder since u(L) is renormalized 
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to zero. In the remaining region we have to distin- 
guish the cases K u ^ K w . For K w < K < K u (i.e. 
qp < 4a/3) and uo = the phase correlations are lost 
on the scale of the KT correlation length £ w of the vor- 
tex unbinding transition. Not too close to this transition 
£ W A « e (Scarc/2h)(i-K v /K) _ Switching on the disorder, 
u will be renormalized by strong phase fluctuations which 
lead to an exponential decay of u ~ Uoe - constL /t™ such 
that disorder is irrelevant for the vortex gas as long as 
£iu £u ■ We expect that the relation £ w sa £ u determines 
the position of the phase boundary between a pinned low- 
K phase where vortices are favoured by the disorder and 
an unpinned high-i-T phase where vortices are induced 
by quantum fluctuations. This line terminates in K w for 
Mo — > (see Fig. 2). If Score is large, £ w will be large 
as well and £ w w £ u will be reached only for K ss K*. 
For moderate values of Score, the unpinning transition 
may be lowered considerably by quantum phase slips. In 
the opposite case K u < K < K w (i.e. qp > 4y3) phase 
fluctuations renormalize weak disorder to zero such that 
vortices are still suppressed until K reaches K w where 
vortex unbinding occurs. In this case two sharp phase 
transitions have to be expected. 

Our flow equations describe also the effect of a com- 
mensurate lattice potential on the CDW: if the wave 
length ir/kp of the CDW modulation is commensurate 
with the period a of the underlying lattice such that 
ir/kp = n/(qa) with n, q integer, an Umklapp term 
w cos qtp appears in the Hamiltonian Q . We obtain 
the results in this case from (^) - (^|) (and the con- 
clusions derived from them) if we use the replacements 
K -»■ K~ l ,t -> t/K 2 and w -> w/K 2 . Thus the lattice 
potential is relevant for K < K w with — 8/q 2 . 

Next we consider the application of the results ob- 
tained so far to a one- dimensional Bose fluid. Its density 
operator is given by eq. (1) if we identify Q/ir = po = p\ 
(p = 2). d x ip is conjugate to the phase 9 of the Bose field 
|||. With the replacements K -> K' 1 , t -> t/K 2 and 
w = 0, (^) describes the action of the lD-supcrfluid in 
a random potential, v denotes the phase velocity of the 
sound waves with v/(irK) — pa/m and irvK = k/(tt 2 Pq) 
where k is the compressibility. The transition between 
the superfluid and the localized phase occurs for K* = 
2/3 H . Thermal fluctuations again suppress the disorder 
and destroy the superfluid localization transition in ID. 
In contrast to CDWs here the (9-field may have vortex- 
like singularities in space-time and the flow equations 
a PPly again. The vortex unbinding transition ap- 
pears at K w with JQ, = 1/2, (q = 2). If both w and 
u are non-zero we can use the canonical transformation 
to rewrite the vortex contribution in the form w cos(qip). 
For K < K U ,K W both perturbations are irrelevant and 
the system is superfluid. For K w < K < K u the de- 
cay of u is stopped due to the suppression of the if fluc- 



tuations due to w. An Imry-Ma-argument shows fur- 
ther, that the q(p = ir(n + 1/2) state is destroyed on the 
scale £ w £ w /u 2 (log{;) by arbitrary weak disorder, i.e. 
vortices become irrelevant above this scale. On larger 
scales one can expects that quantum fluctuations wash 
out the disorder, the system is still superfluid. Finally, 
at K > K U ,K W both perturbations are relevant and su- 
perfluidity is destroyed. 

To conclude we have shown, that in ID CDWs/SDWs 
and superfluids disorder driven zero temperature phase 
transitions are destroyed by thermal fluctuations leav- 
ing behind a rich cross-over behavior. Quantum phase 
slips in CDWs and superfluids lead to additional phase 
transitions and shift the unpinning transition in CDWs 
to smaller if-values. Coulomb hardening and dissipative 
quantum effects will be discussed in a forthcoming pub- 
lication hj. 
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